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We present an analytic treatment near the phase transition for the critical temperature of (3+1)- 
dimensional holographic superconductors in Einstein-Gauss-Bonnet gravity with backreaction. We 
find that the backreaction makes the critical temperature of the superconductor decrease and con- 
densation harder. This is consistent with previous numerical results. 
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I. INTRODUCTION 

Gauge/Gravity duality stems from string theory provides a rich tool for analyzing strongly coupled field theory pj. 
Especially, the duality provides an established method for calculating correlation functions in a strongly interacting 
field theory using a dual classical gravity description 0, Q . Holographic superconductors established in [H, H| are 
remarkable examples where the Gauge/Gravity duality plays an important role. There, a superconducting phase 
transition is described by black hole physics according to the duality. 

Many of previous works on the holographic superconductors were performed in the probe limit, where the backre- 
action of matter fields on the spacetime metric is neglected. However, we found that the effect of the Gauss-Bonnet 
coupling lowers the critical temperature of holographic superconductors in previous work Q . Backreaction becomes 
important when considering lower temperature of black holes in AdS spacetime, because lower Hawking temperature 
of black holes means smaller black holes, i.e., larger Coulomb energy of the matter fields near the black hole horizon. 

The critical temperature was obtained numerically both with and without the backreaction [1, 0, • As an analytic 
approach for deriving the critical temperature, an approximate analytic formula was proposed in the probe limit by a 
matching method [6(. An alternative analytic method using the expansion around the critical point where the phase 
transition occurs was also proposed Q . Some other analytic approaches were also proposed in the probe limit [9l-fl2| . 
In this note, we derive analytically the critical temperature of the Gauss-Bonnet holographic superconductors with 
backreaction by combining the small backreaction approximation and the matching method near the critical point. 

II. GAUSS-BONNET BLACK HOLES 

We begin with the action for a Maxwell field and a charged complex scalar field coupled to the Einstein-Gauss- 
Bonnet: 
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where g is the determinant of a metric g^ v and R pv \ p , R pv and R are the Ricmann curvature tensor, Ricci tensor, 
and the Ricci scalar, respectively, q is the charge of the scalar field. We take the Gauss-Bonnet coupling constant 
a to be positive. Here, the negative cosmological constant term —6/L 2 is also introduced. We look for electrically 
charged plane-symmetric hairy black hole solutions taking the metric ansatz 



dr 



ds 2 = -f(r)e 2u{r) dt 2 + ^ + ^{dx 2 + dy 2 + dz 2 ) 



f(r) L 2 

together with a static ansatz for the fields, 

A» = (#r), 0,0, 0,0), 



ip = ip(r). 
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Here, without loss of generality ip can be taken to be real. First, we look for the solutions in normal phase, ip = 0. 
Wc find v is constant and 
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Here, (i and p arc interpreted as a chemical potential and charge density of the dual theory on the boundary, 
respectively. We also find 
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where we chose the minus sign of the solutions so that we have a solution in the Einstein limit (a — > 0). The horizon 
radius r + is defined through the requirement that f{r+) = 0, so we set a constant of integration by imposing this 
condition at the horizon. We see the above solution becomes the Reissner-Nordstrom-AdS solution in the Einstein 
limit. In order to avoid a naked singularity, we need to restrict the parameter range as a < L 2 / '4. For general a, the 
solution (O behaves as 
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in the asymptotic region. Hence, we define the effective asymptotic AdS scale by 
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Here, the maximum value of a is called the Chern-Simons limit. The Hawking temperature is given by 
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where a prime denotes derivative with respect to r. This will be interpreted as the temperature of the holographic 
superconductors. 



III. GAUSS-BONNET HOLOGRAPHIC SUPERCONDUCTORS 



In order to obtain the solutions in superconducting phase, i/i ^ 0, we need to take into account the boundary 
conditions at the horizon and the AdS boundary. The position of the horizon, r+, is defined through /(r+) = 0. We 
can set v(r+) = there by rescaling the time coordinate. Then the Einstein equations give: 



( r+ )2 e 2„(r + ) 



Regularity at the horizon for </> and ifj gives two conditions: 



As we want the spacetime to be asymptotically AdS, we look for a solution with 
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0(r+) = 0, = y^-y^(r + ) . (11) 



u(r) = const. , f(r) = ^ , (12) 

at the AdS boundary. 

• Asymptotically (r — > oo) the solutions of and ^ are found to be: 

^~^- + ^, (13) 

where A± = 2 ± ^/4 + m 2 L 2 . Note that these are not entirely free parameters, as there is a scaling degree of freedom 
in the equations of motion. As in we impose that p is fixed, which determines the scale of the system. For tp, in 
order to have a normalizablc solution we take C_ = 0. 
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According to the Gauge/Gravity duality PHH, we can interpret (Oa + ) = C+, where Oa + is the operator with 
the conformal dimension A + dual to the scalar field. Thus, we are going to calculate the condensate (Ca+) for fixed 
charge density. We note that since C_ is regarded as the source term of the operator, we put it zero at the end by 
using the Gauge/Gravity dictionary @, Q, which is consistent with taking C_ = 0. 

Let us change the coordinate and set z — r + jr. Under this transformation, the Einstein, Maxwell and the scalar 
equations become 
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where the prime now denotes a derivative with respect to z. The region r + < r < oo now corresponds to < z < 1. 
If one sets </> = <f>/q, tp = ip/l m t ne action (fTJ), the Maxwell and the scalar equations remain unchanged, while the 
gravitational coupling of the Einstein equations changes k 2 — > k 2 /q 2 . If one takes the limit q — > oo, the matter sources 
drop out of the Einstein equations and this is the probe limit. To go beyond the probe limit, we can take either finite 
q with setting 2k 2 = 1, or finite 2k 2 with setting q = 1. The paper [j| took the former choice to consider the effects 
of backreaction of the spacetime metric, but we will take the latter choice in the following. 

In order to solve these equations, we focus on near the critical point as in 0, [l3[, around which the stability was 
confirmed (l4l |. It is convenient to introduce a scalar operator as an expansion parameter: 



(18) 



As -0 is small near the critical point, we expand tp from the first order. From Eq. (|16p . 4> an d ip arc expanded by e 2 
subsequently as follows: 
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and in this situation where starting from the normal phase, the background can be expanded around the Reissner- 
Nordstom-AdS spacetime: 



f = fo + e 2 / 2 + e 4 /4 + • • ■ , 
v = e 2 i>2 + e 4 ^4 ■ ■ • , 



(21) 
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where £<1. 

While ji is also expanded by the order parameter as 



fx = fj, + e 5^2 
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where §fj,2 is also positive. So, we find the order parameter as a function of the chemical potential, which is expressed 
by 
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The exponent 1/2 is consistent with the Ginzburg-Landau mean field theory for phase transitions. When fj. = [io, the 
order parameter becomes zero, which means the critical value of fj, is defined by /Lto such as \i c = /io- Now we begin 
to solve equations order by order. 

At zeroth order, we impose 4>o(l) = at the horizon, z = 1, and 4>o(0) = Ho = const, at the boundary, z = 0. Then 
the solution of Eq. (fl"6)) is given by 



<p (z) = /i (l - z 2 ) 
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This gives a relation p = /io r + by the coordinate transformation. Inserting this solution into Eq. (|15[) . we obtain 
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where we chose the minus sign of the solutions so that we have a solution in the Einstein limit. We also used /o(l) = 
at the horizon as we did in Eq. . We find the above solution becomes the Rcissncr-Nordstrom-AdS solution in the 
Einstein limit. 

Even at first order, it is difficult to deal with both of the backreaction on the spacetime metric and the Gauss-Bonnet 
term. Here we will appeal to the matching method used in [(|. Before proceeding to it, let us check the regularity 
condition at the horizon at this order: 
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The behavior of tJj at the boundary is given by 
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where A± is the same as in Eq. (|13l) . For the boundary conditions, we take D_ to be zero as we did in Eq. (flU)) . 

Now we find the solution of tpi under these conditions by the matching method and derive the critical temperature. 
We can expand ipi in a Taylor series near the horizon as: 
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Here, we have Eq. (|2T|) for the first order coefficients of i/ji, and without loss of generality we can take V'i(l) > to 
have ip\{z) positive. The second order coefficients of tp± are computed by using Eq. (|17l) . Then we find the second 
derivative at the horizon is expressed by 
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After eliminating tp[(l) from above equation by using Eq. (|27j). an approximate solution near the horizon is given by 
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Here, V'i(l) is still unknown. 

On the other hand, from (|28[) . ipi in the asymptotic region are given by 



ipx(z) ~ D + z 
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We have set D_ = from the boundary condition. Here D + is another unknown constant. 

Now we try to connect the solutions Eq. (|3"Tj) and (f3"2"j) smoothly at z m in order to obtain V'i(l) and D+. As we shall 
see below the choice of z m does not change the qualitative features of solutions. In order to connect those solutions 
smoothly, we require the following two conditions: 



A + z^- 1 D + =' + 



rim 



Vi(l)-2 



4/6(1) 



4/6(1 
, , fo (1) 



./6(1) ^(l); 1 4 ./6(1) 2 



(l-^)Vi(l), (33) 



./6(i) /6(i) J 4 /6(i) 



(i-£ ro )V>i(i) 



From Eqs (|33|) and (f34|) . we find the relation between ipi (1) and £)_| 
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Plugging this back into Eq. (f34|) . we find the following relation in order to get a non-trivial solution, V'i(l) 7^ 0, 
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Similarly, plugging Eq. (f35|) back in Eq. (f33|) we can get the solution of ipi (1) and then D + as well in principle. But 
what we want to know here is the critical temperature rather than deriving the solution of tjj up to this order, so we 
are going to focus on it in the following. For this purpose, putting the values of fb(l), fo (1) and 0q(1) in the above 
relation, then Eq. Q36p yields the equation for /jq: 
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In order to solve the above equation with respect to /io, we assume k 2 < 1 in the following. This means that all 
functions are expanded by k 2 as well. That is, solutions near the phase transition are obtained in the small backreaction 
approximation together with the matching method. As the first term disappears from the above equation, we can 
solve it easily and we get 
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where /^o is positive. Combining Eq. (|38[) with the relation //o = p/ r +i which is given under Eq. (|25p. we find r + is 
given by 
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The Hawking temperature Eq. (JSJ) up to this order becomes 

/^(iK«(D _ r+ 



32 

5 — 4z m m 2 L 2 



m L — (1 — z rn )m a 
4m 2 



-1/6 



1 - Z„ 



18 



-Q 



Iff = -' 



47rra 



7tL 2 



1-2^1 



(39) 



(40) 



The critical temperature is defined at the point where the order parameter becomes zero, which leads to Tji = T c at 
Ho = He- Eliminating Ho{ = He) from the above temperature Tjj using Eq. (|38[) . we get the critical temperature 
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Substituting r + with Eq. 



, we obtain the critical temperature of the form 
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(41) 
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IV. CONCLUSION 



We presented an analytic treatment near the phase transition for the critical temperature of (3 + l)-dimcnsional 
holographic superconductors in Einstein-Gauss-Bonnct gravity with backreaction by matching the solution expanded 
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from infinity with that expanded from the horizon. In this method, there is an ambiguity in the choice of the matching 
radius. However, the result turns out to be fairly insensitive to the choice of it. The result reproduces our previous 
analytic calculation for m? = — 3/L 2 , k 2 = with the choice z m — 1/2, which was in good agreement with the 
numerical result @. It also agrees with 0] for m 2 = —A/L 2 , k 2 = a = with the choice z m = 1/2. It may be noted 
that the choice z m ~ 0.5 is roughly equal to z = \/r+/L, corresponding to the geometrical mean of the horizon radius 
and the AdS scale, r = yJr + L. We found that the coefficient of the corrections due to backreaction, T 2 , is positive 
definite irrespective of the value of a. This means the effects of the backreaction makes condensation harder. This 
result also agrees with the numerical results obtained in (l5l - fl7j . 
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